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ABSTRACT 

LIN,  L.;  DEMIRBILEK,  Z.,  and  MASE,  H.,  2011.  Recent  Capabilities  of  CMS-Wave:  A  Coastal  Wave  Model  for 
Inlets  and  Navigation  Projects.  In:  Roberts,  T.M.,  Rosati,  J.D.,  and  Wang,  P.  (eds.),  Proceedings,  Symposium  to 
Honor  Dr.  Nicholas  C.  Kraus,  Journal  of  Coastal  Research,  Special  Issue,  No.  59,  pp.  7-14.  West  Palm  Beach 
(Florida),  ISSN  0749-0208. 


The  Coastal  Inlets  Research  Program  (CIRP)  of  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC), 
Coastal  and  Hydraulics  Laboratory  (CHL)  has  developed  a  nearshore  spectral  wave  transformation  numerical  model 
to  address  needs  of  the  U.S.  Army  Corps  of  Engineers  (USACE)  navigation  projects.  The  model  is  called  CMS- 
Wave  and  is  part  of  Coastal  Modeling  System  (CMS)  for  wave  estimates  in  the  vicinity  of  coastal  and  estuarine 
navigation  channels.  It  can  simulate  important  wave  processes  at  coastal  inlets  including  wave  diffraction,  refraction, 
reflection,  wave  breaking  and  dissipation  mechanisms,  wave-current  interaction,  and  wave  generation  and  growth. 
This  paper  describes  recent  improvements  in  CMS-Wave  that  include  semi-empirical  estimates  of  wave  run-up  and 
overtopping,  nonlinear  wave-wave  interactions,  and  wave  dissipation  over  muddy  bottoms.  CMS-Wave  may  be  used 
with  nested  grids  and  variable  rectangular  cells  in  a  rapid  mode  to  assimilate  full-plane  wave  generation  for 
circulation  and  sediment  transport  models.  A  brief  description  of  these  recent  capabilities  is  provided. 
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INTRODUCTION 

Improved  and  practical  wave  models  are  needed  for  USACE 
navigation  mission  dealing  with  wave  transformation  in  the 
vicinity  of  coastal  inlets,  estuaries  and  river  mouths,  harbors, 
and  navigation  channels.  Wave  interaction  with  coastal  jetties, 
breakwaters  and  revetments  requires  modeling  wave  reflection, 
diffraction,  transmission  and  overtopping  for  accurate  wave 
estimates  of  channel  infilling  and  potential  breaching  or  flanking 
of  structures.  Calculations  are  required  for  wave  propagation 
across  channels  as  well  as  over  or  through  structures,  where 
combined  wave  diffraction  and  transmission  are  present 
simultaneously.  Wave  run-up  and  overtopping  of  structures  are 
frequently  needed.  To  address  these  needs,  a  spectral  wave 
transformation  model  called  CMS-Wave  has  been  developed  to 
support  the  operation  and  maintenance  of  coastal  inlet  and 
navigation  projects,  also  the  risk  and  reliability  assessment  of 
shipping  in  inlets  and  harbors.  It  is  a  steady-state  half-plane 
wave  transformation  model,  where  wave  energy  propagation  is 
independent  of  time  and  calculated  from  seaward  towards  the 
coastal  boundary.  Wave  reflection  off  coastal  structures  and 
beaches  is  included  in  CMS-Wave.  The  workings  and 
performance  of  model  have  been  described  by  Lin  et  al.  (2006 
and  2008),  Seabergh  et  al.  (2008),  and  Demirbilek  et  al.  (2009). 
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The  focus  of  this  paper  is  to  provide  a  brief  description  of  some 
recent  capabilities  added  to  the  model. 

CMS-Wave  can  represent  important  coastal  wave  processes 
including  diffraction,  refraction,  reflection,  wave  breaking  and 
dissipation  mechanisms,  wave-current  interaction,  and  wave 
generation  and  growth.  The  wave  diffraction  is  calculated  in  the 
wave-action  balance  equation  by  an  additional  term.  Recent 
features  incorporated  for  grid  nesting,  variable  rectangle  cells, 
wave  run-up  and  overtopping,  and  assimilation  of  full-plane 
wave  generation  are  to  support  circulation  and  sediment 
transport  studies.  Details  of  the  theory  and  numerical 
implementation  are  presented  in  the  CMS-Wave  technical  report 
(Lin  et  al .,  2008),  including  a  number  of  examples  of  practical 
applications.  In  this  paper,  we  provide  validation  for  wave  run¬ 
up,  transmission  and  overtopping  of  structures  with  laboratory 
data,  wave  dissipation  over  muddy  beds  and  nonlinear  wave- 
wave  interaction  with  the  field  and  experiment  data. 

WAVE-ACTION  BALANCE  EQUATION 

CMS-Wave  calculates  the  spectral  wave  transformation  based 
on  the  wave-action  balance  equation  (Lin  et  al .,  2008): 

dCN  dC  N  dCnN 
— ^  +  — 1— + — §— = 

dx  dy  dQ 

( CCg  cos2  QNy  )^  -  cos2  QNyy  +Sln+Sdp+Snl 

(1) 
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where  N  =  E  /  o  is  the  frequency  and  direction  dependent  wave- 
action  density,  defined  as  the  wave  energy-density  E  =  E  (x,  y, 
o,  0)  divided  by  the  intrinsic  frequency  a.  Nv  and  Nyy  denote  the 
first  and  second  derivatives  with  respect  to  y;  x  and  y  are  the 
horizontal  coordinates;  0  is  the  wave  direction  measured 
counterclockwise  from  the  x-axis;  C  and  Cg  are  wave  celerity 
and  group  velocity;  Cx ,  Cy,  and  C0  are  the  characteristic  velocity 
with  respect  to  x,  y,  and  0,  respectively,  k  is  an  empirical 
parameter  representing  the  intensity  of  wave  diffraction  effect. 
The  right-hand  side  terms  respectively  are:  Sin  is  the  source  (e.g., 
wind  input),  Sdp  is  the  sink  (e.g.,  bottom  friction,  wave  breaking, 
whitecapping,  etc.),  and  Sn/  is  the  nonlinear  wave -wave 
interaction. 

The  first  term  on  the  right  side  of  Equation  1  is  the  wave 
diffraction  term  formulated  from  a  parabolic  approximation 
wave  theory  (Mase,  2001).  In  applications,  the  diffraction 
intensity  parameter  k  (>  0)  needs  to  be  calibrated  and  optimized 
for  featured  structures.  The  model  omits  the  diffraction  effect  for 
k  —  0  and  calculates  the  diffraction  for  k  >  0.  In  practice,  the 
value  of  k  may  range  from  0  (no  diffraction)  to  4  (strong 
diffraction)  for  calculating  diffraction  effects.  A  constant  value 
of  k  —  2.5  has  been  used  by  Mase  et  al.  (2001,  2005a,  2005b)  to 
simulate  wave  diffraction  for  narrow  and  wide  gap  breakwater 
applications.  Lin  et  al.  (2008)  and  Demirbilek  et  al.  (2009) 
demonstrate  that  value  of  k  =  4  is  appropriate  for  semi-inifmite 
long  breakwaters  and  also  in  narrow  gaps  (inlets)  with  openings 
equal  or  less  than  one  wavelength.  For  wider  gaps  with  the 
opening  greater  than  one  wavelength,  k  =  3  is  recommended. 
The  exact  value  of  k  in  an  application  is  dependent  on  the 
structure’s  geometry,  local  bathymetry  and  incident  wave 
conditions,  and  may  need  to  be  fine-tuned  with  data.  The  default 
value  of  k=  4  is  used  in  the  model,  corresponding  to  strong 
diffraction.  Implementation  of  wave  diffraction  is  approximate, 
and  phase-resolving  wave  models  (Holthuijsen  et  al.,  2004)  may 
be  used  to  verify  estimate  of  waves  near  structures,  inlets  and 
harbors. 

Figures  1  and  2  show  an  example  of  the  CMS-Wave 
computational  domain  and  calculated  wave  results  at  four 
different  locations  (HMB1  to  HMB4),  respectively,  at  Half- 
Moon  Bay  in  Grays  Harbor  entrance,  Washington,  USA. 
Osborne  and  Davies  (2004)  described  the  field  data  collection  in 
Grays  Harbor.  Numerical  simulation  is  conducted  for  December 
10  through  31,  2003.  The  input  wind  and  wave  data  were 
obtained  from  the  National  Data  Buoy  Center  (NDBC  - 
http://www.ndbc.noaa.gov)  Station  46029  and  Coastal  Data 
Information  Program  (http://cdip.ucsd.edu)  Buoy  036  (NDBC 
46211),  respectively.  CMS-Wave  was  run  in  the  coupled  mode 
with  a  CMS-Flow  model  (Buttolph  et  al.,  2006; 
http://cirp.usace.army.mil/products)  that  provided  calculated 
water  level  and  current  input  fields.  The  effect  of  current  on 
waves  (Hs  shown  in  Figure  2  is  the  significant  height  defined  as 
the  mean  of  the  highest  1/3  wave  height  )  is  pronounced  at 
gauges  HMB1  and  HMB2  located  closer  to  the  navigation 
channel  in  relatively  deep  water.  Comparison  of  simulation 
results  with  the  field  data  shows  that  shallow  water  effects  on 
wave  diffraction,  refraction,  and  breaking  are  evident  at  HMB3 
and  HMB4. 


Figure  1.  CMS-Wave  grid  and  data-collection  stations  at  Grays  Harbor, 
WA.  USA. 


Figure  2.  Measured  and  calculated  waves  at  HMB1  to  4,  10-31 
December  2003. 
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WAVE  RUN-UP 

The  wave  run-up  consists  of  two  components:  (a)  a  rise  of  the 
mean  water  level  resulting  from  wave  breaking  at  the  shore, 
known  as  the  wave  setup  q,  and  (b)  swash  oscillation  of  the 
incident  waves.  The  wave  setup  in  CMS-Wave  is  computed 
from  the  horizontal  momentum  equations  by  neglecting  current, 
surface  wind  drag  and  bottom  stresses  as 


while  Mase  and  Iwagaki  (1984)  specified  the  incident  wave  by 
JONSWAP  spectrum  (Hasselmann  et  al .,  1973)  as  follows: 

Ej{a)  = 

^-j-exp  |-1 .25(-^-)4 1  exp  jin  y  exp[-0.5(^^=-)2  ]  j 

(9) 


&n=  1 

I  (2) 

8x  pgh 

^  dx 

fy  j 

1 

Sn=A_  I 

(3) 

8y  pgh  \ 

v  dx 

J 

where  p,  g ,  h  are  the  water  density,  gravitational  acceleration 
and  water  depth,  respectively.  Three  radiation  stress  components 
Sxx,  Sxy ,  Syy  are  calculated  using  the  linear  wave  theory  (Dean 
and  Dalrymple,  1984)  as 


s«  =  J 

\E 

n(cos20+l)--^ 

<76  do 

(4) 

\E 

n(sin20  +  l)-| 

<70  do 

(5) 

S„=j 

|  ^-n  sin  20  d6  do 

(6) 

where  _±  kh  and  k  is  the  wave  number.  The  swash 

2  sinh2A:/z 

oscillation  of  waves  on  the  beach  face  is  a  random  process.  In 
engineering  applications,  the  2-percent  runup  defined  as 
exceedance  of  the  vertical  level  is  commonly  used  and  denoted 
as  R2%  (or  R2).  Komar  (1998)  gives  an  estimate  of  R2%  on 
beaches,  seawalls  and  jetties  as 

R2%=2\i!\  (7) 


The  R2%  is  calculated  using  Equation  7  at  the  land-water 
interface  and  averaged  with  the  local  depth  to  determine  if  water 
can  flood  the  dry  cells.  If  the  wave  run-up  level  is  higher  than 
elevation  of  adjacent  dry  land  cells,  these  cells  are  flooded  to 
simulate  the  overtopping  and  overwash  at  these  cells. 

Wave  run-up  calculation  in  CMS-Wave  is  verified  with  data 
collected  from  two  laboratory  studies  of  random  wave  up-rush 
on  plane  smooth  slopes  (Ahrens  and  Titus,  1981;  Mase  and 
Iwagaki,  1984).  There  are  together  a  total  of  395  run-up  cases 
from  these  two  studies  covering  a  broad  range  of  plane  slope 
(1:1  to  1:30)  and  incident  spectrum.  Ahrens  and  Titus  (1981) 
characterized  the  incident  wave  by  the  spectral  peak  frequency 
and  a  wave  group  parameter  (Goda,  1970)  defined  as 


where 


e„  = 


<jE2(<j)d<j 

[\E(a)dcjf 


E(cj)  =  |  E(cr ,  6)d  6 


(8) 


where  a ,  p  are  scale  and  shape  parameters,  and  crm  ,  y  are 

spectral  peak  frequency  and  peak-enhancement  factor, 
respectively.  Incident  waves  in  Ahrens  and  Titus  are  converted 
to  a  JONSWAP  spectrum  as  input  to  CMS-Wave  using  the 
relationship  ofy  =  2(Q  -3) +  3  derived  from  Equations  8  and  9. 

Figure  3  shows  the  calculated  and  measured  R2%.  Overall, 
calculated  run-up  estimates  agree  with  data;  the  correlation 
coefficient  for  all  cases  is  0.83  (Lin  et  al. ,  2008).  For  steeper 
slopes  (1:1  to  2:3),  the  model  tends  to  overpredict  R2o/o 
suggesting  Equation  7  is  more  accurate  for  mild  slopes. 


Figure  3.  Measured  and  calculated  2%  exceedence  wave  runup. 


WAVE  TRANSMISSION  AND  OVERTOPPING  OF 
STRUCTURES 

Wave  transmission  over  low-crested  or  submerged 
impermeable  structures  is  exhibited  by  fall  of  the  overtopping 
water  mass.  The  ratio  of  structure  crest  elevation  to  incident 
wave  height  is  a  key  parameter  governing  the  wave 
transmission.  CMS-Wave  calculates  the  transmission  coefficient 
Kh  defined  as  the  transmitted  wave  height  divided  by  the 
incident  wave  height,  based  on  the  simple  expression  (Goda, 
2000): 
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K.  =0.3  (1.5 — — ),  for  0<  — <1.25  (10) 

H,  H, 

where  h  is  the  crest  elevation  of  the  breakwater  above  the  still- 

c 

water  level  ( hc  is  negative  for  a  submerged  breakwater),  and  Ht 

is  the  incident  wave  height.  For  a  composite  breakwater, 
protected  by  a  mound  of  armor  units  at  its  front,  Kt  is  calculated 
(Goda,  2000)  as 

K.  =0.3  (1.1— L)  for  0<  — <0.75  (U) 

t  TT  J '  TT 


Wave  overtopping  rate  is  calculated  using  two  different 
methods.  In  a  coupled  CMS -Wave  and  CMS -Flow  simulation 
( e.g .,  flow  model  gets  radiation  stresses  from  wave  model  and 
provides  calculated  current  and  water  levels  to  wave  model),  the 
overtopping  rate  is  the  flow  rate  calculated  by  CMS-Flow. 
Without  coupling,  CMS-Wave  calculates  the  average 
overtopping  rate  using  the  empirical  formula  of  Hughes  (2008) 
given  by 


q  =  JgHf 


0.0336  +  0.53 


(13) 


For  permeable  rubble-mound  breakwaters,  the  transmission  is 
calculated  using  d’Angremond  et  al.  (1996)  formula  as 


K,  =  0.64 


1  -  exp  (A)]  (^)-°'31  -  0.4^,  for  B  <  10  Ht 

(12) 


where  B  is  the  crest  width  and  £  is  the  Iribarren  parameter 

defined  as  the  fore-slope  of  the  breakwater  divided  by  the 
square-root  of  deepwater  incident  wave  steepness.  In  practice, 
Equations  10  to  12  are  applicable  to  both  monochromatic  and 
random  waves.  Figure  4  shows  the  calculated  transmission 
coefficients  and  data  curves  compiled  by  Goda  (2000)  for  a 
vertical  breakwater.  In  these  simulations,  CMS-Wave  was 
forced  by  a  monochromatic  wave  of  1  m  and  6  sec  for  a  vertical 
breakwater  with  /z  =  10  m,  d  =  5  m,  and  B  =  20  m.  Table  1 
shows  comparison  of  calculated  transmission  coefficients  and 
laboratory  data  compiled  by  Goda  (2000),  and  Equations  10  and 
1 1  for  monochromatic  incident  waves. 


Figure  4.  Measured  and  calculated  transmission  coefficients  for  a 
vertical  breakwater. 


Table  2  provides  comparison  of  the  calculated  average 
overtopping  rates  obtained  with  decoupled  and  coupled  CMS 
runs  performed  for  laboratory  measurements  of  Hughes  (2008) 
for  four  incident  irregular  waves  and  two  steady  storm  surges 
overtopping  a  levee.  In  the  simulation  of  higher  surge  (water 
level  =  1.3  m)  and  larger  waves  (incident  wave  height  =  2.3  m), 
the  calculated  overtopping  rates  from  the  CMS-Wave  alone 
agree  better  with  data  than  results  from  the  coupled  run.  This 
suggests  that  empirical  formulas  used  in  CMS-Wave  for 
overtopping,  runup  and  transmission  need  further  verification 
with  data  to  determine  appropriateness  of  each  formula. 


Table  1.  Comparison  of  transmission  coefficients,  Kt 


}ic  (m)  Vertical  Breakwater  Composite  Breakwater 


Data  (Goda,  2000)  Eg,  10  CMS-Wave  Eg,  11  CMS-Wave 


-2 

1.02 

1.02 

-1.5 

1.03 

1.03 

-1 

0.76 

0.78 

0.78 

-0.5 

0.60 

0.63 

0.63 

0 

0.42 

0.45 

0.46 

0.33 

0.34 

0.5 

0.26 

0.30 

0.27 

0.18 

0.18 

1 

0.13 

0.15 

0.15 

0.03 

0.04 

1.5 

0.06 

0.10 

0.024 

2 

0.05 

0.07 

0.018 

Table  2.  Comparison  of  average  overtopping  rate  (m2/sec) 


Exp  No.  Surge 
Level 
(m) 

Wave 

Height 

(m) 

Wave 

Period 

(sec) 

Measured 
(Hughes,  2008) 

CMS 

Calculated 

CMS-Wave 

R128 

0.29 

- 

0.27 

0.28 

- 

0.82 

6.1 

0.38 

0.38 

0.39 

R109 

0.29 

- 

- 

0.26 

0.28 

- 

2.48 

13.7 

0.70 

0.85 

0.92 

R121 

1.3 

- 

- 

2.55 

2.57 

- 

2.3 

6.1 

2.67 

2.93 

2.76 

R127 

1.3 

- 

- 

2.54 

2.57 

- 

2.3 

14.4 

2.84 

2.98 

2.81 

WAVE-WAVE  INTERACTIONS 

The  exact  solution  of  wave-wave  interaction  requires  solving 
computationally  expensive  six-dimensional  integral  (Resio  and 
Tracy,  1982).  The  wave-wave  interactions  in  CMS-Wave  are 
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calculated  using  a  modified  method  of  Jenkins  and  Phillips 
(2001).  They  proposed  a  simple  formula  to  represent  wave-wave 
interactions  as  a  second-order  diffusion  operator  of  the  isotropic 
form  that  conserves  the  wave  action.  CMS- Wave  implements 
this  formula  directly  in  the  wave-action  balance  equation  such 
that  no  additional  integration  is  required  to  calculate  wave-wave 
interactions.  By  keep  only  the  first  and  second-order  terms  in  the 
Jenkins-Phillips  formulation,  wave-wave  interactions  in  finite 
water  depth  can  be  expressed  as 


=  a 


dF  ,  d2F 

- +  b — - 

da  dO 


(14) 


where  *  =  A[i+(2«-l)2cosh2tt]-l  is  a  function  of  kh,  b=j_, 
2  n  ncj 


and 


F  —  k3a5 


(2  nfg 


(15) 


The  variation  of  a  and  bcr  as  functions  of  kh  is  shown  in 
Figure  5,  indicating  that  wave -wave  interactions  are  more 
significant  in  the  intermediate  water  depths  and  diminish  in  the 
shallow  water.  This  is  because  wave-wave  interactions  initially 
take  place  in  the  deepwater,  continue  to  evolve  in  the 
intermediate  depth,  and  gradually  diminish  in  shallow  water. 
The  effect  of  interactions  on  wave  evolution  is  greater  over  long 
fetches  in  a  large  ocean  domain  and  less  in  a  local  coastal 
region.  Figure  6  shows  the  comparison  of  directionally 
integrated  Snl  from  Equation  14  and  exact  computations  in  the 
examples  of  Hasselmann  et  al.  (1985).  These  calculations  are  for 
the  JONSWAP  spectra  with  y  =2  and  y  =5.  Figure  7  shows  the 

calculated  nonlinear  wave  energy  transfer  rate  Snl  in  the 
frequency  and  direction  domain  for  a  JONSWAP  spectrum  with 
y  =5.  These  calculated  results  are  consistent  with  the  observed 

and  theoretical  results  that  the  nonlinear  wave -wave  interactions 
cause  wave  energy  to  transfer  from  high  to  low  frequencies 
(downshifting). 

Figure  8  shows  the  calculated  wave  height  fields  with  and 
without  nonlinear  wave-wave  interactions  at  the  Louisiana  coast 
in  the  north  central  Gulf  of  Mexico.  The  input  wind  and  incident 
wave  data  are  supplied  by  NDBC  Buoy  42041  located 
approximately  in  the  middle  of  offshore  boundary.  Figure  9 
shows  the  corresponding  wave  period  fields  with  and  without 
the  nonlinear  wave  energy  transfer.  The  simulation  is  for  the 
growth  of  an  incident  deepwater  southeast  wave  of  0.3  m  and  4 
sec  under  a  moderate  15-kt  wind  from  southeast.  With  the 
nonlinear  wave -wave  interactions,  the  calculated  wave  height 
along  the  coast  is  about  10  percent  higher  than  without  the 
nonlinear  energy  transfer  effect.  The  nonlinear  wave-wave 
interactions  also  increase  wave  periods  over  a  larger  coverage 
area,  which  is  a  consequence  of  wave  energy  transfer  from 
higher  to  lower  frequencies.  Wind  input  was  triggered  in  these 
simulations.  It  should  be  noted  that  wind  input  can  reveal  better 
wave -wave  interaction  result  because  more  wave  energy  from 
wind  input  can  transfer  from  high  to  low  frequencies.  The  wave- 


Figure  5.  Nonlinear  wave-wave  interaction  coefficients  a  and  ba  as 
functions  of  kh. 


Figure  6.  Comparison  of  directionally  integrated  S„i  for  JONSWAP 
spectrum  with  y  =2  and  5. 


wave  interaction  can  be  rather  small  and  insignificant  in 
magnitude  without  the  wind  input. 

The  wave  transformations  near  inlet  jetties  and  in  the  wave 
diffraction  zone  can  be  improved  with  the  nonlinear  wave -wave 
interactions.  Figure  2  shows  the  comparison  of  calculated  waves 
with  and  without  the  nonlinear  energy  transfer  in  Grays  Harbor. 
The  calculated  wave  height  with  the  nonlinear  energy  transfer  is 
about  5  percent  higher  than  without  the  nonlinear  energy 
transfer.  This  change  is  small  and  in  practical  applications  can 
be  neglected.  More  importantly,  we  note  that  the  proposed  new 
formulation  that  extends  the  Jenkins-Phillips  original 
formulation  from  deepwater  to  finite  depth  should  be 
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Figure  7.  Calculated  non-linear  energy  transfer  for  a  directional 
JONS  WAP  spectrum  with  y  =5. 


Figure  8.  Calculated  wave  height  field  (a)  with  and  (b)  without  non¬ 
linear  wave  energy  transfer. 


extensively  tested  in  different  depths  with  different  wave  and 
wind  conditions.  The  robustness  and  consistency  of  our 
proposed  formula  must  be  validated  with  laboratory  and  field 
data  and  confirmed  against  other  formulations. 

WAVE  DISSIPATION  OVER  MUDDY  BOTTOMS 

The  wave  dissipation  over  muddy  beds  in  CMS- Wave  is 
calculated  based  on  the  assumption  that  the  turbulent  eddy 
viscosity  is  several  orders  of  magnitude  greater  than  the 
kinematic  viscosity  of  sea  water.  By  neglecting  the  kinematic 


Figure  9.  Calculated  wave  period  field  (a)  with  and  (b)  without  non¬ 
linear  wave  energy  transfer. 


viscous  effect,  the  wave  dissipation  over  a  muddy  bed  can  be 
expressed  as  (Lamb,  1932) 


Sdp=-Myk+vt)k2E  (16) 


where  vk  is  the  kinematic  viscosity  and  vt  is  the  turbulent  eddy 
viscosity.  In  the  present  study,  vt  is  formulated  by  a  maximum 
viscosity  yfb  representing  the  wave  breaking  condition  times  the 

ratio  of  wave  height  over  depth  as  v  _  v  Hs_ .  The  value  of  Vtb 

1  tb  h 

is  set  to  0.04  m2/sec  for  the  muddy  bed  and  0.01  m2/sec  for  the 
primary  sand  area  based  on  model  comparison  with  data.  As  an 
example,  wave  fields  affected  by  muddy  beds  for  the  Louisiana 
coast  are  calculated  and  shown  in  Figures  8  and  9.  In  Figure  10, 
CSI3  is  the  location  of  muddy  bed,  and  CSI5  is  for  the  primary 
sandy  bed.  These  wave  gauges  were  deployed  by  the  Coastal 
Studies  Institute  (CSI)  at  Louisiana  State  University  (Sheremet 
and  Stone,  2003;  Sheremet  et  al.,  2005).  Figure  11  shows  the 
comparison  of  calculated  and  measured  wave  heights  at  gauges 
CSI3  and  CSI5.  The  calculated  wave  heights  at  the  muddy  bed 
CSI3  using  Equation  16  agree  better  with  measurements  as 
compared  to  those  without  the  muddy  bed  effect  included.  For 
the  more  sandy  location  CSI5  away  from  the  muddy  area,  the 
calculated  wave  heights  are  less  affected  by  inclusion  of  the 
muddy  bed  effect  .The  present  treatment  of  muddy  bed  does 
address  the  potential  affect  of  mud  on  dynamics  of  waves  as 
reported  by  other  researchers  (Sheremet,  20005;  Kaihatu, 
2008).  The  goal  in  the  present  study  is  to  provide 
approximate  estimate  of  mud  on  spectral  wave  dissipation, 
without  using  any  complicated  approaches. 
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Figure  10.  Location  of  muddy  bed,  NDBC  42041,  and  CSI  wave  gauges. 


Figure  11.  Measured  (+)  and  Calculated  ( — )  waves  on  muddy  bed  along 
Louisiana  coast. 


CONCLUSIONS 

This  paper  introduces  some  recent  capabilities  that  have  been 
added  to  CMS-Wave  model.  Due  to  space  limitation,  only  a 
brief  description  of  the  added  features  is  provided  and  details  of 
these  will  be  reported  soon  in  technical  reports  and  notes. 
Examples  provided  for  the  new  capabilities  are  simply 
comparisons  to  analytical  solutions,  other  models  and  available 
data.  The  comparisons  include  wave  run-up,  wave  transmission 
over  and  through  structures,  wave-wave  interactions,  and  wave 
energy  dissipation  on  muddy  beds.  These  new  features 
implemented  in  CMS-Wave  may  not  be  essential,  but  often  are 
needed  for  wave  estimates  in  coastal  engineering  applications 
near  structures  and  in  navigation  channels,  where  wave-current 
and  wave-structure  interactions  can  become  important  and 
significant.  Additional  research  is  in  progress  to  further  evaluate 
the  appropriateness  of  the  proposed  capabilities,  which  require 
extensive  validation  and  verification  before  these  new  features 
can  be  confidently  used  in  the  engineering  practice.  Once 
completed,  a  comprehensive  mathematical  description  of  these 
new  added  features  will  be  provided  with  validation  examples. 


Other  features  can  also  be  important  concerning  wave 
asymmetry  in  shallow  water,  infra-gravity  waves  for  seiching  in 
harbors,  and  surface  roughness  of  breakwaters.  Future 
companion  papers  in  this  series  will  address  a  more  complete 
description  of  these  nearshore  wave  processes. 
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